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Abstract. 

We present a holographic description of four-dimensional single-scalar inflationary universes 
in terms of a three-dimensional quantum field theory (QFT). The holographic description 
correctly reproduces standard inflationary predictions in their regime of applicability. In the 
opposite case, wherein gravity is strongly coupled at early times, we propose a holographic 
j-h ' description in terms of perturbative QFT and present models capable of satisfying the current 

observational constraints while exhibiting a phenomenology distinct from standard inflation. 
This provides a qualitatively new method for generating a nearly scale-invariant spectrum of 
primordial cosmological perturbations. 
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1. Introduction. 

The notion of holography emerged from black hole physics as an answer to the question: why 
is the entropy of a black hole proportional to the area of its horizon? Since entropy is an 
extensive quantity, one would have expected the entropy to be instead proportional to the 
volume the black hole occupies. Typically, entropy is a measure of the number of degrees of 
freedom. The scaling of the gravitational entropy has thus been taken as an indication of a 
new fundamental principle, the holographic principle, that underlies any quantum theory of 
gravity. More precisely, holography states that any quantum gravitational system in (d + 1) 
dimensions should have a dual description in terms of a QFT without gravity in one dimension 
less [1]. Indeed, if gravity is holographic this would explain the scaling behavior of the black 
hole entropy, since the entropy of a QFT scales like the volume, which is the same as the area 
in one dimension higher. 

Holography provides a new paradigm for physical reality, the consequences of which we are 
only just beginning to comprehend. According to this picture, one of the macroscopic dimensions 
of spacetime and one of the forces in the universe, namely gravity, are emergent phenomena in 
an underlying lower-dimensional QFT. Concrete realizations of holography have been found in 
string theory, and a precise holographic dictionary was established shortly thereafter [2, 3, 4]. 
To date, almost all such realizations involve spacetimes with a negative cosmological constant. 
The arguments that led to the holographic principle apply more generally, however, and suggest 
that one should be able to establish a holographic dictionary that applies to our own universe. 
The purpose of the work presented in [5], and further discussed here, is to propose a concrete 
holographic framework that applies to our own universe, and in particular to its cosmological 
evolution. More precisely, we will describe how to set up holography for inflationary cosmology. 

Any holographic proposal for cosmology should specify (i) what the dual QFT is, and (ii) 
how it can be used to compute cosmological observables. Having defined such a duality, the new 
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description should recover established results in their regime of applicability. Indeed, we will see 
that our holographic models correctly reproduce standard inflationary results when standard 
inflation is applicable, namely, under the assumption that gravity was weakly coupled at early 
times. Yet perhaps more importantly, our holographic approach also gives a qualitatively new 
method for generating a nearly scale-invariant spectrum when gravity was strongly coupled at 
early times. As we will discuss in the following, there exist holographic models that are capable 
of satisfying all current observational constraints while exhibiting a phenomenology distinct from 
standard inflation. 

Over the last two decades striking new observations have transformed cosmology from a 
qualitative to quantitative science [6]. These observations reveal a spatially flat universe, 
endowed with small-amplitude primordial perturbations that are approximately Gaussian and 
adiabatic with a nearly scale-invariant spectrum. This data is consistent with the generic 
predictions of inflationary cosmology and set inflation as the leading theoretical paradigm for 
the initial conditions of Big Bang cosmology. Yet inflation, despite these successes, is still 
unsatisfactory in a number of ways: it generically requires fine tuning and there are trans- 
Planckian issues and questions about the initial conditions for inflation [7]. 

With future observations promising an unprecedented era of precision cosmology, it becomes 
imperative that inflation is embedded in a UV complete theory (indeed there is increasing amount 
of effort devoted to embedding inflation in string theory) , and it is also important that alternative 
scenarios are developed. The holographic approach that we undertake provides both. Firstly, 
the holographic models we will discuss here are three-dimensional super-renormalizable theories 
and thus are UV complete. Secondly, holographic dualities are strong/weak coupling dualities 
meaning that in the regime where one description is weakly coupled, the other is strongly 
coupled. This provides an arena for constructing new models with intrinsic strong-coupling 
gravitational dynamics at early times that possess only a weakly coupled three-dimensional 
QFT description, and are thus outside the class of model described by standard inflation. As 
we will see, such models may lead to qualitatively different predictions for the cosmological 
observables that will be measured in the near future. Furthermore, quite apart from the fresh 
perspective on early universe cosmology such an approach offers, there are also a number of 
more pragmatic reasons for developing a holographic framework for cosmology: uncovering the 
structure of three-dimensional QFT in cosmological observables brings in new intuition about 
their structure and may lead to more efficient computational techniques, cf. the computation of 
non-Gaussianities in [8]. 

The holographic description we propose uses the one-to-one correspondence between 
cosmologies and domain-wall spacetimes discussed in [9, 10] and assumes that the standard 
gauge/gravity duality is valid. More precisely, the steps involved are illustrated in Fig. 1. The 
first step is to map any given inflationary model to a domain-wall spacetime. For inflationary 
cosmologies that at late times approach either a de Sitter spacetime or a power-law scaling 
solution 1 , the corresponding domain- wall solutions describe holographic renormalization group 
flows. For these cases there is an operational gauge/gravity duality, namely one has a dual 
description in terms of a three-dimensional QFT. Now, the map between cosmologies and 
domain- walls may equivalently be expressed entirely in terms of QFT variables, and amounts to 
a certain analytic continuation of parameters and momenta. Applying this analytic continuation, 
we obtain the QFT dual of the original cosmological spacetime. 

We shall call the resulting theory a 'pseudo'-QFT because we currently have only an 
operational definition of this theory. Namely, we do the computations in the QFT dual to 
the corresponding domain-wall and then apply the analytic continuation. Perhaps a more 

1 This era should then be followed by a hot big bang cosmology, as in standard discussions. Here we only 
discuss the very early universe, i.e., the times when the primordial cosmological perturbations were generated 
(the inflationary epoch) . 
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Figure 1. The 'pseudo'- 
QFT dual to inflationary 
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defined using the corre- 
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fundamental perspective is to consider the QFT action, with complex parameters and complex 
fields as the fundamental objects, and then to consider the results on different real domains as 
applicable to either domain- walls or cosmologies. Note that the supergravity embedding of the 
domain- wall/cosmology correspondence discussed in [11] works in precisely this way. 

This paper is organized as follows. In the next section we discuss the domain- wall/cosmology 
correspondence. Then, in Section 3, we discuss the cosmological observables that we would like 
to compute holographically, and in Section 4 we present the holographic analysis. In Section 5, 
we discuss the analytic continuation to the pseudo-QFT and in Sections 6 and 7 we discuss the 
new models that are strongly coupled at early times but have a weakly coupled QFT description. 

2. Domain- wall/cosmology correspondence. 

We explain in this section the lefthand vertical line in Fig. 1, namely, the correspondence between 
cosmologies and domain-wall spacetimes. For simplicity, we focus on spatially flat universes 
equipped with a single minimally coupled scalar field, but the results can be extended to more 
general cases (eg., non-flat, multi-scalar, non-canonical kinetic terms, etc). The metric and scalar 
field for the unperturbed solution are given by 

ds 2 = ndz 2 + a 2 (z)df 2 , $ = tp(z), (1) 

where rj = — 1 in the case of cosmology, in which case z is the time coordinate, and r\ = +1 
in the case of domain-wall solutions 2 in which case z is the radial coordinate. We take the 
domain-wall to be Euclidean for later convenience. A Lorentzian domain-wall may be obtained 
by continuing one of the x l coordinates to become the time coordinate [10]. The continuation to 
a Euclidean domain-wall is convenient, however, because the QFT vacuum state implicit in the 
Euclidean formulation maps to the Bunch-Davies vacuum on the cosmology side. Other choices 
of cosmological vacuum require considering the boundary QFT in different states, as may be 
accomplished using the real-time formalism of [12]. 

With the appropriate choice of a(z) and tp(z), the configuration in (1) solves the field equations 
that follow from the action 

S = ^j d 4 x^\g~\[-R+(d^) 2 + 2K 2 Vm, (2) 

where k 2 = 8ttG and we are taking the scalar field <3? to be dimensionless. Note that when 
r\ = — 1, the kinetic terms have the appropriate signs for a mostly plus Lorentzian signature 

2 The name 'domain-wall spacetime' dates back to earlier work featuring solutions of this form that interpolate 
between two stationary points of the scalar field potential, one at z = +oo and another at z — -co. In the present 
context the name is somewhat misleading, however, since we consider only the z > part of the geometry. We 
will nevertheless stick with this terminology as it is standard usage in high-energy physics. 
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metric and when rj = +1, the kinetic terms have the correct sign for Euclidean signature. Had 
we expressed both actions in the same signature metric, they would differ only in the sign of 
the potential. It follows that for every flat FRW solution of a model with potential V there is a 
corresponding domain- wall solution of a model with potential — V [9, 10]. 

For background solutions in which the evolution of the scalar field is (piece-wise) monotonic, 
<p(z) can be inverted to give z(cp) permitting the Hubble rate H = a/ a to be re-expressed as 
H(z) = — (1 / 2)W ((p) , where W(<p) is known as the 'fake superpotentiaP. In this case, the 
complete equations of motion for the background take the simple form 



- a = -\w, <p = W >v , 2r ] n 2 V = {W, ip ) 2 -*-W 2 (A) 



This first-order formalism goes back to the work of [13] (for cosmology), where it was obtained 
by application of the Hamilton-Jacobi method. In [10] this formalism was linked to the notion 
of (fake) (pseudo-) supersymmetry. 

We will now extend the correspondence to encompass linear perturbations around the 
background solution. The linearly perturbed metric takes the general form 

ds 2 = r][l + 2<f>(z,x)]dz 2 + 2a 2 (z)[diis(z,x) + Ui(z, x)]dzdx l + a 2 (z)[<% + h ij (z,x)]dx t dx j , 
$ = (p(z)+6(p(z,x), (4) 

where i/j is transverse. The spatial metric perturbation hij may be decomposed as 

hij(z,x) = -2ifi(z,x)5ij + 2didjx(z,x) + 2d { , i w j )(z,x) +jij(z,x), (5) 

where uii is transverse and 7^(2,2:) is transverse traceless. The metric perturbations may then 
be combined into the gauge-invariant combinations 

C = i> + {H/<p)8<p, (6) 

4> = 4> - (<VM, (7) 

v = v -x~ri(Sip/a 2 ip), (8) 
Vi = Vi - coi. (9) 

Physically, £ represents the curvature perturbation on comoving hypersurfaces and has the useful 
property that it tends to a constant on superhorizon scales 3 . 

The equations of motion for cosmological perturbations have been worked out long ago, see 
[15] and references therein, while the corresponding analysis for domain-walls may be found in 
[16, 17, 18]. In the present case, there are two independent perturbations represented by £ and 
7ij, since the Hamiltonian and momentum constraints are equivalent to 

<^=-4 ^=-4i+4' *=°. ( io ) 

H a z H q z 

where q is the comoving wavevector of the perturbations, and the background quantity e(z) is 
defined as e = —H/H 2 = 2{W^/W) 2 . (In standard inflation e would be the usual slow-roll 
parameter; however, we do not assume slow roll here). From the remaining Einstein equations, 
one then finds the following equations of motion for ( and 7^: 

= £ + (3H + e/e)C - r?a~ VC, 

= jij + SHjij - r]a~ 2 q 2 -fij. (11) 

a Following the end of inflation, the superhorizon value of ( then remains constant until horizon re-entry, 
irrespective of the dynamics of the intervening evolution, provided that no entropy perturbations are produced 
[14,15]. 
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Defining now the analytically continued variables R and q according to 

R 2 = — k 2 , q = —iq , (12) 

it is easy to see that a perturbed cosmological solution written in terms of the variables k and 
q continues to a perturbed Euclidean domain-wall solution expressed in terms of the variables 
R and q. Note that (11) only requires q 2 = —q 2 and in (12) we made a choice of a branch cut in 
the function q = \Jlp (for reasons to be explained in the next section). As it is clear from (3), 
continuing k is equivalent to continuing V. Here we prefer to continue k since, as we will see, 
the former has a clear interpretation in terms of dual QFT variables. 

We have thus established that the correspondence between cosmologies and domain-walls 
holds, not only for the background solutions, but also for linear perturbations around them. 
This is the basis for the relation between power spectra and holographic 2-point functions, to 
be discussed momentarily. The argument can be generalized to arbitrary order to relate non- 
Gaussianities to holographic higher-point functions [19]. 

3. Cosmological observables. 

In the inflationary paradigm, cosmological perturbations originate on sub-horizon scales as 
quantum fluctuations of the vacuum. Quantizing the perturbations in the usual manner, one 
finds the scalar and tensor superhorizon power spectra 

A|W = |^(C(g)C(-g)) = |^IC 9 (o)l 2 , 

A|(<7) = ^<7tf(9)7«(-9)> = ^l7 g( 0)| 2 , (13) 

where J q (o) and Q(o) are the constant late-time values of the cosmological mode functions 7 g (z) 
and ( q {z). 

The mode functions are themselves solutions of the classical equations of motion (11) (setting 
7*i = lq e iji f° r some time-independent polarization tensor e^). To select a unique solution 
for each mode function we impose the Bunch-Davies vacuum condition Cqilq ~ exp(— igr) as 
r — > — oo, where the conformal time r = f z dz'/a(z'). The normalization of each solution (up 
to an overall phase) may then be fixed by imposing the canonical commutation relations for the 
corresponding quantum fields. This leads to the Wronskian conditions, 

i = c,nf )* - <) C *, i/2 = 7(?n w* - nj%*, (14) 

where U { q ° = 2ea 3 K~ 2 C q and nj 7) = (l/4)a 3 k 2 7q are the canonical momenta associated with 
each mode function, and we have set h to unity. 

To make connection with the holographic analysis to follow, we introduce the linear response 
functions E and O satisfying 

n<p = fK 9 , nW = £ 7? . (is) 

These quantities are well-defined since we have already selected a unique solution for each mode 
function. Substituting these definitions into the Wronskian conditions, which are valid at all 
times, the cosmological power spectra may be re-expressed as 

^' tMnW ' A ^ )= 2,qX)M ' <16) 
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where Imf2( ) and ImE^ are the constant late-time values of the imaginary part of the response 
functions (more precisely, the subscript indicates that this is the part of the response function 
that is invariant under dilatations, see the discussion at the end of subsection 4.1). We will see 
shortly how the response functions also give the 2-point function of the pseudo-QFT. 

Let us now consider the corresponding domain-wall solution obtained by the applying the 
continuation (12). The early-time behavior ~ exp(— iqr) of the cosmological perturbations maps 
to the exponentially decaying behavior ~ exp(gr) in the interior of the domain-wall (r — > — oo). 
Such regularity in the interior is a prerequisite for holography, explaining our choice of sign in 
the continuation of q. 

The domain- wall response functions E and Q. [18] are defined analogously to (15), namely 

n< c) = -nc,- u^ = -E 7q , (17) 

where njp = 2eo?R~ 2 Qq and IL^ = {l/4)o?R~ 2A jq are the radial canonical momenta. The minus 
sign in (17) are inserted so that 

n(-iq) = n(q), E(-iq)=E(q). (18) 

By choosing the arbitrary overall phase of the cosmological perturbations appropriately, we 
may ensure that the domain-wall perturbations are everywhere real. The domain-wall response 
functions are then purely real, while their cosmological counterparts are complex. 

4. Holographic analysis. 

In this section we briefly review relevant material from gauge/gravity duality, corresponding to 
the upper horizontal line in Fig. 1. There are two classes of domain- wall solutions for which 
holography is well understood: 

Asymptotically AdS domain-walls. In this case the solution behaves asymptotically as 

a(z) ~ e z , <p ~ as z — > oo. (19) 

The boundary theory has a UV fixed point which corresponds to the bulk AdS critical point. 
Depending on the rate at which <p approaches zero as z — > oo, the QFT is either a deformation 
of the conformal field theory (CFT), or else the CFT in a state in which the dual scalar 
operator acquires a nonvanishing vacuum expectation value (see [20] for details). Under the 
domain- wall/cosmology correspondence, these solutions are mapped to cosmologies that are 
asymptotically de Sitter at late times. 

Asymptotically power-law solutions. In this case the solution behaves asymptotically as 

a(z) ~ (z/zo) n , (p ~ v / 2nlog(z/zo) as z — > oo, (20) 

where zq = n—1. This case has only very recently been understood [21]. For n = 7 the 
asymptotic geometry is the near-horizon limit of a stack of D2 brane solutions. In general, 
these solutions describe QFTs with a dimensionful coupling constant in the regime where the 
dimensionality of the coupling constant drives the dynamics. Under the domain- wall/cosmology 
correspondence, these solutions are mapped to cosmologies that are asymptotically power-law 
at late times. 
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4-1. Basics of holography. 

Gauge/gravity duality is an exact equivalence between a bulk gravitational theory and a 
boundary QFT. Typically, the boundary QFT is a gauge theory that admits a large N expansion. 
The N here denotes the rank of the gauge group: an example of such theory, with gauge group 
SU(N), is discussed in section 7.1. The large N limit consists of taking N — > 00 while keeping 
fixed the 't Hooft coupling constant A = gy M N [22]. One can show that in this limit only planar 
diagrams survive. On the bulk side, taking the large N limit means that one suppresses loop 
effects. The value of A then controls whether the supergravity approximation is valid or not. 

The duality relates bulk fields with gauge-invariant operators of the boundary theory. For 
example, the bulk metric corresponds to the boundary stress-energy tensor, Tj, while bulk 
scalar fields correspond to boundary scalar operators, such as tr i^-i™, where Fij is the field 
strength of the gauge field and the trace is over the gauge group indices. Correlation functions 
of such gauge invariant operators can then be extracted from the asymptotics of bulk solutions 
and conversely, given correlation functions of the dual operators, one can reconstruct asymptotic 
solutions. 

To understand the holographic computations we need to know a few things about the 
structure of asymptotic solutions of the field equations. The discuss below refers specifically 
to a four dimensional bulk spacetime (the general features are the same in any dimension but 
the details are different). The results for the case of asymptotically power-law spacetimes can 
be obtained from the results for asymptotically AalS2 a +i spacetimes via a dimensional reduction 
on a T 2<T_3 torus and analytic continuation in a [23]. The most general asymptotic solution for 
both cases can be shown to be of the form [24, 21] 

ds 2 = dr 2 + gij(r, x)dxMx J , 
gij (r, x) = e 2r (g (0)ij (x) + e- 2r g {2)ij (x) + • • • + e~ 2(7r g {2a)ij (x) + • • •) , (21) 

where a = 3/2 for solutions of Einstein's equations with negative cosmological constant and 
o~ = (3n — l)/2(n — 1) > 3/2 for the case of asymptotically power-law solutions. In the latter 
case, the metric above is given in the so-called dual frame [25]: the Einstein frame metric 
gfj in (2) is related to the metric g® in (21) via the Weyl transformation gf- = exp(\&)g® , 

where A = \/2/n. In (21), gto)ij(. x ) is an arbitrary (non-degenerate) three-dimensional metric, 
from which the g(2k)ij( x ) with k < a are locally determined, while g^2a)ij( x ) i s onr y partially 
constrained by the asymptotic analysis of the field equations. One can show, however, that 
this coefficient is directly related to the expectation value of the boundary stress-energy tensor 
[24,21]: 

(Tij) = ^ ( 2 ^(2^). (22) 

Consider now a scalar field 4 ^ that is dual to an operator O of dimension A. The conformal 
dimension depends on the mass of ^ via m 2 = A (3 — A) in the asymptotically AdS case and is 
equal to A = 4 in the asymptotically power-law case (for appropriately normalized O, see [21]). 
The asymptotic expansion for ^ has a form analogous to (21), 

r ) = e ( A - 3 >(vI/ (0) + e - 2 ^ (2) + • • • + e- 2 ^ (25) (x) + •••), (23) 

Here VP (q) 0*0 is unconstrained and is the source for the dual operator 0, with all subleading 
terms, up to order exp(— 2ar), then being locally determined in terms of the sources. The 
following term in the series, *&(2a)i is undetermined and is related to the expectation value of 

4 In the asymptotically AdS case, = <f>, while in the asymptotically power-law case, = exp((n — l)$/\/2n). 
To match with the discussion in [21] for the D2-brane case (n = 7), note that $hero = — (\/l4/5)</>there, where 
^there IS the scalar field in [21] for p = 2. 



The Holographic Universe 



8 



dual operator (O). In the asymptotically AdS case, 2a = (2A — 3), while in the asymptotically 
power-law case, a = a = (3n — l)/2(n — 1). (For details, see [20, 21]). 
The constraints on g(2a)ij due to the Einstein equations imply 

V l (T lJ ) + (O)V.^ {0) =0, (24) 
<2?) + (3-A)* (0 )<O)=0, 

and these are precisely the expected Ward identities 5 . 

The relation (21)-(22) can be read in two ways: (i) given a supergravity solution it allows us to 
read off the QFT data encoded by the solution; (ii) given QFT data it provides a reconstruction 
of the bulk spacetime in the neighborhood of the boundary. Note that the latter is true even 
when the supergravity approximation is not valid in the interior of spacetime (because the 
curvature is large there). The terms exhibited in (21), apart from g^a)ijj are non-normalizable 
terms and are not affected by dynamics. The first term that is affected by dynamics is g(2o)ij 
and this is indeed unconstrained asymptotically, except for the constraints due to symmetries 
(namely the Ward identities (24)). The gauge/gravity duality provides a dual description of 
dynamics, so the statement of the duality is that gi2 a )ij determined from QFT via (22) should 
agree with the value obtained from string dynamics in a spacetime with these asymptotics. When 
the gravity approximation is valid throughout, the asymptotics yield sufficient information to 
uniquely reconstruct a regular bulk solution: the pair (<7(o)jj, 9(2a)ij) (and similar for the scalar 
field) are coordinates in the covariant phase space of the gravitational theory [26]. 

An alternative way to express these results that we will use in the next subsection is to use 
the radial Hamiltonian formulation of [18]. This is a Hamiltonian formulation where the radial 
direction plays the role of time. In this formalism the expectation value of the stress-energy 
tensor is given by 

(T ij ) = (^U ij ) , (25) 
\V9 /(3) 

where fly is the radial canonical momentum and we will momentarily explain the meaning of 
subscript (3). There is a similar formula for the expectation value of O that we will not need 
here. A fundamental property of the spacetimes (21) is that the radial derivative is, to leading 
order as r — > oo, proportional to the dilatation operator So, 

d r = 8 D (l + 0[e~ 2r }), (26) 

where 

5 D 9ij{x,r) = 2 gij (x,r), 5 D ^(x,r) = (A - 3)tf(z,r) (27) 

(one can verify (26) by inspection of (21) and (23)). The expressions in (27) are the standard 
dilatation transformation rules for a metric and a source that couples to an operator of dimension 
A in three spacetime dimensions. Equation (26) is a precise version of the often quoted relation 
between the radial direction and the energy scale of the dual QFT. In our context, equation (26) 
implies that one can trade the radial expansion for an expansion in terms of eigenfunctions of 
dilatation operator. An eigenfunction of weight n is by definition, 

5 D A (n) = -nA (n) . (28) 

As follows from (26), Ar n \ ~ e _nr (l + 0[e _2r ]), so the radial expansion and the expansion is 
eigenfunctions of the dilatation operator are closely related. The chief advantage of the radial 

5 Note there is no conformal anomaly in three dimensions. 
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Hamiltonian reformulation is that the expansion is now manifestly covariant, whereas expanding 
in a particular coordinate is not a covariant operation. Now, the radial canonical momentum 
can be decomposed in eigenfunctions of the dilatation operator and the subscript (3) in (25) 
indicates that we should pick the part with dilatation eigenvalue 6 3. One may have expected 
this on general grounds, since the dimension of stress energy tensor is three dimensions is 3. 
More generally, the expectation value of an operator of dimension A is given by the piece of 
the corresponding radial canonical momentum of weight A. Note that in general the radial 
canonical momentum contains also pieces with weight less than A (i.e., less than 3 in the case 
of Tj): these pieces diverge as r — > oo. One of the advantages of radial Hamiltonian formalism 
[18] is that renormalization is simple: one simply discards all pieces with weight less than A 
(plus additional logarithmic divergences at weight A, when present). One can indeed show that 
removing these pieces is equivalent to adding local covariant counterterms to the on-shell action. 
The radial canonical momentum for the asymptotically AdS case is equal to 

U ij = ^^g(K ij -Kg ij ), (29) 

with Kij = {l/2)d r gij the extrinsic curvature of constant-r slices. In the case of asymptotically 
power-law backgrounds, the relevant canonical momentum is that of the dual frame [21], namely 

= Vse A * {Kij — (K + \d r $)6*) , (30) 

where A = y/2/n and the metric and extrinsic curvature are those of the dual frame. 



4-2. 2-point functions. 

To linear order in the sources, the variation of the 1-point function for the stress-energy tensor 
is 

S(Tij(x)) = - j d 3 Vy ^ (^(T tJ (x)T kl (y))5g^(y) + {Tij(x)0{y))5 m (y)) . (31) 

It follows that in order to obtain the holographic 2-point functions we need to solve the 
the linearized equations of motion about the domain-wall solution with Dirichlet boundary 
conditions at infinity and imposing regularity in the interior. Transforming to momentum space, 
on general grounds, the 2-point function takes the form 

{Tij{q)T kl {-q)) = A(q)Uijkl + B(q)ir ijm , (32) 

where Hijki is the three-dimensional transverse traceless projection operator defined by 

n ijfcZ = T.{^ik^lj + nuTT kj - 1Tijir k i), -Kij = 5ij - ^If . (33) 

z q 
Decomposing the boundary metric as in (5), we then find 

5{Tij) = \A{q) lm - 2B{q)^ {0) T H j - (Tij(q)0(-q))8 m . (34) 

From (25), this expression is to be compared with the bulk radial canonical momentum expanded 
to linear order. We begin by writing the perturbed metric as in (4) and (5), with the lapse and 

6 In odd (bulk) dimensions, the transformation rule of this specific coefficient has also an additional anomalous 
contribution due to the conformal anomaly. There is no conformal anomaly in our case, and this coefficient is a 
true eigenfunction of the dilatation operator. 
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shift perturbations gauged to zero (<f> = v = Vi = 0), and we additionally set u>i to zero using 
the constraints (10) and a spatial gauge transformation. Then, for the case of asymptotically 
AdS backgrounds, using (29) we find 



5{Ti) = ^[5K5) 



^](3) 



1 



2^8) + fx*) + ^7 



1 

2^ 



(3) 



Expressing (10) in the present gauge, it follows that 

2ip = ip8(p, 



(35) 
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where in the last line we have used the definition of the domain-wall response function Q in (17) 
and expanded £ according to (6). Then, using the definition of the response function E also in 
(17), we find 
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Since the scale factor a has dilatation weight —1 (as follows from (27)), comparison with (34) 
yields 7 

1 



A(q) = 4E (0) (q) 



B(q) = -n (0) (q), 



(38) 



where the zero subscript indicates the pieces of the response functions that have zero weight 
under dilatations, so in particular they are independent of r as r — > oo. Note that, in general, E 
and diverge as r — > oo. Extracting £7 ) and fj( ) correctly then requires first determining the 
terms with eigenvalue less than zero and subtracting these from E and f2, before taking the limit 
r — > oo (see [18] and the example in the next subsection). The issue here is that the subtraction 
of the infinite pieces may induce a change in the finite part as well. This can happen if the local 
covariant counterterms needed to cancel the infinities necessarily have a finite part as well. 

Considering now the case of backgrounds that are asymptotically power-law, as before 
we decompose the metric perturbations as in (4) and (5), and gauge the shift and vector 
perturbations to zero (u = = oji = 0). This time, however, we choose the lapse perturbation 
to be 4> = (\/2)8ip. Transforming to the dual frame, we then find 



ds 2 = e-^ds 2 = dr 2 + a 2 [bij + h^dx'dx 3 , 
hij = -2^8ij + 2didjX + 7ij> 



(39) 



where a = exp(— X(p/2)a, with A=y2/n and the radial coordinate dr = exp(— \(p/2)dz. The 
dual frame perturbations are related to their Einstein frame counterparts by 



V = V + (X/2)S<p, 
From (25) and (30), we have 



lij — lij- 



(40) 



5(7]) = ^ [e^ ({8K + \8v, r )8) - 8K) + (. . .)6<p 



1 



(3) 



e^ ( 2^ r 8] - \8<p, r 8] + fx,r^j + \l\r + (• • -)6<P 



J (3) 



s 3\tp/2 



2^8)+q 2 X *) + \i l j + {---)5v 



J (3) 



(41) 



7 The term proportional to q 2 in the coefficient of ip in (37) only contributes a contact term to B(q) which we 
drop. 
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It follows from (34) that the terms proportional to Sip do not contribute to the stress tensor 
2-point function that we are interested in, so in the above (and below) these terms have been 
suppressed. The Einstein frame constraints (10), when expressed in the present gauge, yield 



2^ = (. . .)Sip, 



q 2 x 



+ 



R 2 n 



a?H 2a 3 



^ + (...)<*¥>■ 



Thus we have 



6(Ti) 



(42) 



(43) 



J (3) 



and since the dilatation weight of a is —1, we again recover (38) modulo contact terms. 



4-3. Example: Power-law inflation. 

To illustrate the above discussion, let us consider the domain-wall backgrounds equal (rather 
than asymptotic) to (20) discussed in [21], namely 



a = (z/z ) n , ip = V2nln(z/z ), 



n 



1 > 0. 



(44) 



Under the domain- wall/cosmology correspondence these solutions are mapped to cosmologies 
undergoing exact power-law inflation. While this model is strongly constrained by the WMAP 
data [6], this need not concern us here since our purpose is simply to illustrate the steps involved 
in the holographic computation. Furthermore, in Section 7 we will see that the strong coupling 
version of these models (i.e., where gravity is strongly coupled at early times but the dual 
three-dimensional QFT is weakly coupled) are compatible with observations. 

Following [10], one can obtain the fake superpotential from the solution yielding W = 
— (2n/zo)ex.p(—ip/V2n). It follows that e = 1/n and both 7^ and £ obey the same equation of 
motion, which for the domain-wall spacetime reads 



= C g -+(3n/z)C 9 --(z/z )" 2 VC 
Imposing regularity in the interior, the solution is 



(45) 



(46) 



where K c is a modified Bessel function of the second kind of order a = (3n — l)/2(n — 1) > 3/2, 
Cq is an arbitrary function of q and the radial coordinate p = q{z / Zo) l ~ n . The boundary z — > 00 
corresponds to p = while the domain- wall interior corresponds to p — > 00. The corresponding 
radial canonical momentum is equal to 



n. 



(0 



2e 



Expanding about p = 0, we find 



Cq = Cn 1 + 



1 



2Cg ( P 



-2(7 



pd p (p*K a (p)). 



(47) 



9 r(l-cr) 9 „ 



n 



(0 



q nR 2 \2{l - a)' 4<T(1 + a) 



and thus 



n(q) = - 



n 



(0 



2q 



2(7 



(q nR 2 \2(l-a) 



.2(1-.) _ 2ar(l - a) 
J 4-r(l + a) 



(48) 



(49) 
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As expected, this diverges as p — > 0. To compute the 2-point function we need to identify the 
parts that have negative dilatation eigenvalue, subtract them from (49) and then take p — > 0. 

To do this, we first transform to the dual frame via gfj = e 

-y/{2/n)v g E and then c hange radial 
variable, r = zqIti^z/zq) = —\n(p/q). The metric is now that of AdS, 

ds 2 = dr 2 + e 2r dx 2 , (50) 

and the dilatation operator is exactly equal to the radial derivative, 

5 D = d r = -pd p . (51) 

(This reflects the fact that the AdS isometry group is the same as the conformal group in one 
dimension less). It follows that any monomial in p is an eigenfunction of Sd, 

5 DP n = -np n , (52) 

and one can simply identify in (49) all terms with negative eigenvalue; for example, ^(_2 CT +2) = 
g 2<T /(nK 2 (l-a))p- 2 ( <T - 1 ). We then have 

In this example, the identification of the terms with negative eigenvalues could be 
accomplished by inspection. In more complicated examples, however, this is no longer the case, 
so we briefly indicate here how one could compute them (see [18] for a more complete discussion). 
Starting from (45), and using the definition of the corresponding canonical momentum (the first 
equality in (47)), one obtains 

I%- - ^< = 0. (54) 

Inserting in this equation the definition of the response function, and changing the radial 
coordinate from z to r, one obtains 

d r n - a 2 £l 2 e- 2 ° r + ( f/a^e 2 ^- 1 ^ = 0, (55) 

where a 2 = nR 2 /2. This may now be solved asymptotically by expanding f2 in dilatation 
eigenvalues, 

^ = ^2&(-2*+2k), (56) 

k>l 

making use 8 of d r = Sd and collecting all terms with the same weight. For example, to leading 
order, at weight (— 2<r+2), only the first and last term in (55) can have this weight, and one 
obtains ^(_2 CT +2) = (^ 2 /a 2 2(l— c)) exp(2(<r— l)r) in agreement with (53). Through iteration, 
one may obtain all coefficients with negative eigenvalue. 
Having obtained f2( ), we may finally compute B(q): 

B(q) = ln (0) = - ^ r R 2 Q 2<J = -—wr4 ^ 2 i°- (57) 

w 4 (0) n 4^r(l + cr) y 4 ff r 2 (o-)n sin 7rt H v ' 

8 In examples where the background solution is only asymptotically AdS, the relation between the dilatation 
operator and the radial derivative contains subleading terms (see (26)) and these that must be taken into account 
in this computation, see [18] for a complete discussion. 
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A near-identical argument holds for the tensors 7^ yielding fy ) = (8/n)E< \, and hence 
A(q) = 2nB(q). Via the domain- wall/cosmology correspondence, applying the continuations 
(12) to (53), the imaginary parts of the cosmological response functions are 

4:77 

Imn ( o) = (8/n)Im£ (0) = --^- )K ~ 2 q 2 ° . (58) 
From (16), we then recover the expected cosmological power spectra: 

A|(?) = ^\{q) = n4g "J (a) (59) 

Note that we could equally well have obtained (58) by applying the continuations (12) to 
the unrenormalized domain- wall response function (49), and then taking the imaginary part 
followed by the limit z — > 00. This is because the divergent terms one subtracts to obtain the 
renormalized reponse functions are all analytic functions of q 2 (as may be seen from (49), where 
the leading term is proportional to q 2 ) and hence under the continuation q 2 = —q 2 , these terms 
remain real and do not contribute to the imaginary part of the cosmological response functions. 
Only the leading non-analytic piece of the domain-wall response functions contributes to the 
late-time imaginary part of the cosmological response functions: this leading non-analytic piece 
is finite and is simply Qr \ . In fact, the late-time values of the imaginary parts of the cosmological 
response functions have to be finite as a consequence the Wronskian relations (14) and the fact 
that £ and 7^ tend to a finite constant at late times. 



5. Continuation to the pseudo-QFT. 

We now wish to re-express the bulk analytic continuation (12) in terms of QFT variables, 
corresponding to the vertical line on the right-hand side of Fig. 1. First, the dimensionful 
coupling of the theory, corresponding to the deformation by the operator O can be read off from 
the asymptotics of <3>. Since <£ does not continue, neither does the coupling constant. Second, 
R~ 2 is proportional to the square of the number of colors, N 2 . It follows that the continuation 
(12) amounts to 

N 2 = -N 2 , q=-iq, (60) 

where the barred quantities are associated with the QFT dual to the domain-wall and the 
unbarred quantities are associated with the pseudo-QFT dual to the cosmology. (Thus N is the 
rank of the gauge group of the QFT dual to the domain-wall spacetime, while N is the rank of 
gauge group of the pseudo-QFT dual to the cosmological spacetime). We therefore find that the 
power spectrum for any inflationary cosmology that is asymptotically de Sitter or asymptotically 
power-law can be directly computed from the 2-point function of a three-dimensional QFT via 
the formulae: 

A|(9) = Iftr^Hg) ' AUq) = aiming) " (61) 
This is one of our principal results. 



6. Beyond the weak gravitational description. 

In the discussion so far we have assumed that the description in terms of gravity coupled to a 
scalar field is valid at early times, and that the perturbative quantization of fluctuations can be 
justified. The holographic description also allows us to obtain results when these assumptions 
do not hold. At early times, the theory may be strongly coupled with no useful description 
in terms of low-energy fields (such as the metric and the scalar field). The holographic set-up 
allows us to extract the late-time behavior of the system, which can be expressed in terms of 
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Weakly 

coupled 

gravity 



Strongly 

coupled 

gravity 




Figure 2. Using holography it is 
possible to describe the generation 
of primordial cosmological pertur- 
bations during an early time epoch 
in which the gravitational descrip- 
tion is strongly coupled. At later 
times we envision a smooth tran- 
sition to a conventional hot big 
bang description in which gravity 
is weakly coupled. (Figure adapted 
from [271). 



low-energy fields, from QFT correlators. This is the counterpart of the discussion in Section 
4 where we saw that in gauge/gravity duality the asymptotic behavior of bulk fields near the 
boundary of spacetime is reconstructed by the correlators of the dual QFT. This late-time 
behavior is precisely the information we need to compute the primordial power spectra and 
other cosmological observables. We will assume that this postulated early-time phase in which 
gravity is strongly coupled is subsequently followed by a smooth transition to the usual hot big 
bang cosmology in which gravity is weakly coupled, see Fig. 2. 



7. Holographic phenomenology for cosmology. 

Ideally one would deduce from a string/M-theoretic construction what the dual QFT is. Instead 
we initiate here a holographic phenomenological approach. The dual QFT would involve scalars, 
fermions and gauge fields and it should admit a large N limit. The question is then whether 
one can find a theory which is compatible with current observations. In particular, one might 
consider either deformations of CFTs or theories with a single dimensionful parameter in the 
regime where the dimensionality of the coupling constant drives the dynamics, as these QFTs 
have already featured in our discussion above. 



7.1. A prototype dual QFT. 

We will discuss here super-renormalizable theories that contain one dimensionful coupling 
constant. A prototype example is three-dimensional SU(N) Yang-Mills theory 9 coupled to 
a number of scalars and fermions, all transforming in the adjoint of SU(N). Theories of this 
type are typical in AdS/CFT where they appear as the worldvolume theories of D-branes. A 
general such model that admits a large N limit is 



S 



1 



9ym 



d^xtr 



1 



1 



F^F 1 ^ + ^{D$ J f + ^{D X K f + i> L U)^ L 



(62) 



where we consider Ma gauge fields A 1 (I = 1, . . . , Ma)', Mj, minimal scalars 4> J ( J = 1, . . . , M^); 
M x conformal scalars \ K {K = 1> ■•■) M" x ) and M^ fermions ip (L = 1, M^). The 

couplings XM1M2M3M4 and Mmxi£ (where a and (3 are spinor indices) are dimensionless, and 
we have grouped the scalars appearing in the interaction terms as $ = ({(f) J }, {x K })- Note 



We write the rank of the QFT gauge group here as iV since we will first be performing calculations using the 
QFT dual to the domain-wall spacetime before analytically continuing to the pseudo-QFT. 
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O Figure 3. 1-loop contribution 
to (Tij(q)T k i(-q)}. We sum over 
~'"kl^~ c ^ the contributions from gauge fields, 

scalars and fermions, with each 
diagram yielding a contribution of 
order ^ N 2 q 3 . 



that all terms in the Lagrangian have dimension 4 so this QFT is indeed of the same type as the 
QFTs dual to asymptotically power law solutions. The operator O that featured in our earlier 
discussion of holography is closely related to the Lagrangian (see the discussion in Section 4 of 
[21]). 

The conformally coupled scalars have an R\ 2 coupling when we couple the theory to gravity; 
on a flat background this means the conformally coupled scalars have a different stress energy 
tensor from their minimally coupled counterparts. Specifically, on a flat background, the stress 
tensor is given by 



Tij — 2 tr 
9ym 



2Ft k Fj fc + DtfDrf + D tX K D jX « - -DiD^Y + ^ 7(i ^ 

+ A Mi m 2 m 3 m^ M ^ M ^ M ^ M4 + A lL * M ^ 2 )] ■ (63) 



To extract predictions, we need to compute the coefficients A(q) and B(q) appearing in 
the general decomposition (32) of the stress tensor 2-point function, analytically continue the 
results, and then insert them in the holographic formulae (61) for the power spectra. To 
facilitate comparison with the observational data, we will additionally make use of the standard 
cosmological parameterizations 

/ \ ns(q)— 1 / \ «t(?) 

A 2 s (q) = A 2 s (q ) i^-J , A* (q) = A 2 T (q ) (J- J , (64) 

where A|^ r (go) is the scalar/tensor amplitude at some chosen pivot scale qo, and tis/t(q) is the 
scalar /tensor spectral tilt. 

7.2. 1-loop calculation. 

The leading contribution to the 2-point function of the stress tensor is at one loop (see Fig. 3). 
Since the stress tensor has dimension three, and the only dimensionful quantity that can appear 
to this order is q (1-loop amplitudes are independent of ^yj^f), it follows that 

A(q) = C A N 2 q 3 + 0{g\ M ), B{q) = C B N 2 f + 0( 5ym ), (65) 

where Ca and Cb are numerical coefficients whose value depends only on the field content. 
Explicit calculation then reveals that 

C A = (M A +M 4> +M x + 2A^)/256, C B = (M A + A/^/256. (66) 

Inserting (65) into our holographic formulae (61), we find 

A2 s(l) = m J N 2 CB + O(Sym), A 2 T (q) = + 0(g 2 M ). (67) 
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o 

Figure 4. Diagram topologies contributing at 2-loop order. Each diagram consists of an 
overall factor of N 3 g YM multiplying an integral with superficial degree of divergence two. After 
dimensional regularization and renormalization, the integrals evaluate to ~ Q 2 ^(q/q*), and 
so overall each diagram yields a contribution to the stress tensor 2-point function of order 
~ ^ 3 5YM<? 2 ln(<?/<7*), or equivalents ~ N 2 q 3 g 2 s ki(q/ g*). 




Comparing with (64), we immediately see that the power spectra are scale-invariant to leading 
order (i.e. n$ = 1 + 0(g\ u ), ut = 0(<7 Y m))> regardless of the precise field content of the model. 
To estimate the value of N we may compare with the observed amplitude of the scalar power 
spectrum. From the WMAP data [6] we have A|(c/ ) ~ O(10 -9 ), hence N ~ O(10 4 ), justifying 
our use of the large N limit. 

The observational data also serve to provide an upper bound on the ratio of tensor to scalar 
power spectra. From (67), we find 

r = A^/A| = 32C B /C A , (68) 

and hence an upper bound on r translates into a constraint on the field content of the dual QFT. 
A smaller upper bound on r requires increasing the number of conformal scalars and massless 
fermions and/or decreasing the number of gauge fields and minimal scalars. 



7.3. 2-loop corrections. 

Corrections to the stress tensor 2-point function at 2-loop order 10 give rise to small deviations 
from scale invariance. The full result will be reported elsewhere [19], however, it is easy to obtain 
an order of magnitude estimate on general grounds. The perturbative expansion depends on the 
effective dimensionless coupling constant g 2 s = gy M N/q. Either from inspection or from direct 
calculation of some of the diagrams contributing at 0(g 2 s ) (see Fig. 4), one finds 

A(q) = C A N 2 q 3 [l + D A g 2 cS ln(g/&) + 0(g 4 cS )}, 

B{q) = C B N 2 f[l + D B g 2 cS ln(q/q,) + 0(<&)], (69) 

where Da and Db are numerical coefficients of order one whose value depends only on the field 
content. To compute D A and Db precisely requires summing all the relevant 2-loop diagrams. 
Inserting these two-loop corrected results into the holographic formulae (61), we find 

^ s{q) = le^ivW 1 " Db9 ° s Hq/q * ] + 0{9 ^ 

A2Aq) = " Da9 ^ Hq/q * ] + °^ )] ' (70) 

10 Super-renormalizable theories have infrared divergences, but large N resummation leads to well-defined 
expressions with <?ym effectively playing the role of an infrared regulator. The exact amplitudes are nonanalytic 
functions of the coupling constant [28]. Note that our analytic continuation to pseudo-QFT does not involve the 
coupling constant. 
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Figure 5. The straight line is 
the leading order prediction of 
holographic models with a single di- 
mensionful coupling constant for the 
correlation of the running a s and the 
scalar tilt n s . The data show the 68% 
and 95% CL constraints (marginaliz- 
ing over tensors) at go = 0.002 Mpc -1 , 
and are taken from Fig. 4 of [6]. As 
new data appear the allowed region 
should shrink to a point, which is 
predicted to lie close to the line. 



where the analytically continued effective coupling g 2 s = gy M N/q. In comparison, expanding 
(64) yields 

A 2 s (q) = A|(g )[l + (ns(a)-l) m(g/g ) + O ((n s (g)-l) 2 )], 

A 2 T (q) = A 2 T (q )[l + n T (q) ln(g/g ) + O (n T (g) 2 )]. (71) 

Identifying the renormalization scale g* with the pivot scale go, we then see that the spectral 
amplitudes given in (67) are correct to 0(g^ s ), and that the corresponding spectral tilts are 

n s (q)-l = -D B g 2 s + 0( 5 4 ff ), n T (q) = -D A g 2 s + 0( 5c 4 ff ). (72) 

Comparing with the WMAP data, from Table 4 of [6] we find that (n s -l) ~ O(10~ 2 ) at 
q = 0.002 Mpc -1 , and hence g 2 s ~ O(10~ 2 ) also, justifying our perturbative treatment of the 
QFT. 

To determine whether the spectral tilts are red or blue requires evaluating the signs of Da 
and Db, which will in general depend on the field content of the QFT. It is nonetheless still 
possible to extract predictions which are independent of the field content: for example, in these 
models, the scalar spectral index runs as 

a s = dn s /dlng = -(n s -l) + 0( 5c 4 ff ). (73) 

This prediction is qualitatively different from slow-roll inflation, for which a s /{n s — 1) is of first- 
order in slow-roll [29], yet is nonetheless consistent with the WMAP observational constraints 
on n s and a s given in [6] for a wide range of values of n s and a s , as illustrated in Fig. 5. 

7.4- Non-Gaussianities. 

Once N, g Y M an< ^ ^ ne field content are fixed, all other cosmological observables (such as 
non-Gaussianities, etc.) follow uniquely from straightforward computations. We will present 
details of the correspondence between higher-order QFT correlation functions and non-Gaussian 
cosmological observables elsewhere [19]. Our results indicate, however, that the non-Gaussianity 
parameter /]^| al [30] is independent of iV to leading order, consistent with current observational 
evidence [6]. 
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8. Conclusions. 

Let us summarize the main results. We have presented a holographic framework for early- 
universe cosmology describing the period of time corresponding to the inflationary epoch. In 
particular, we have shown how to compute cosmological observables by performing calculations 
with a three-dimensional dual QFT. This procedure was discussed explicitly for the case of 
the primordial power spectrum, which is related to the 2-point functions of the dual QFT. 
Higher-point functions are related to non-Gaussianities, as will be discussed elsewhere. When 
gravity is weakly coupled at early times, holography correctly reproduces standard inflationary 
predictions for cosmological observables. When gravity is instead strongly coupled at early 
times, one finds new models that have a weakly coupled QFT description. We saw how models 
of this type exist that are compatible with current observations yet nevertheless have a distinct 
phenomenology from standard inflation. The proposed holographic approach thus provides 
a qualitatively new method for generating a nearly scale-invariant spectrum of primordial 
cosmological perturbations. 

A special case to consider is de Sitter spacetime. A correspondence between de Sitter (dS) and 
CFT has been proposed in the past [31], and it is natural to wish to understand the relation of 
the present work to dS/CFT. Under the domain- wall/cosmology correspondence, dS is mapped 
to AdS, and our results for dS follow via suitable analytic continuation from the results for 
AdS. It is well known that quantum correlators, such as, for example, the Feynman propagators 
of massive scalar fields in AdS and dS spacetime, do not map to each other under analytic 
continuation, see for example [32]. This has been one of the obstacles in trying to establish 
a dS/CFT correspondence using analytic continuation. In our case, however, we do not map 
quantum de Sitter correlators to quantum AdS correlators. Instead, we map dS correlators 
directly to correlation functions of the dual QFT, and (as we have seen) one can successfully 
establish such a map. An approach to dS/CFT that has more in common with the present 
work is that of Maldacena [8], although our analytic continuation is different: we analytically 
continue both the momenta and Newton's constant, whereas in [8] it was the dS radius that 
was continued. Furthermore, compared to previous works [33, 34, 35, 36] that focused mostly 
on the computation of the scalar power spectrum for asymptotically dS geometries, we gave a 
complete discussion of both the scalar and tensor power spectra and our work includes the case 
of asymptotically power-law inflation. Finally, we proposed a precise definition of the dual QFT. 

While uncovering an underlying holographic structure in inflationary cosmology is 
conceptually important, this in itself would lead to no new hard results were one to remain 
in the regime where inflationary computations are well-justified. The reason is that in this 
regime, the holographic computations simply reproduce well-known results. Rather, the power 
of the holographic approach is that it leads to new models in which the gravitational dynamics 
were strongly coupled at early times; for these models the standard inflationary computations 
do not apply. This is precisely the distinguishing feature of the present work: because we give 
an explicit definition of the dual QFT, we are able to obtain just such models with strongly 
coupled gravitational dynamics at early times. These models have a weakly coupled dual QFT 
description, permitting the analysis of a scenario that would otherwise be quite intractable. 

We have seen how the dual QFT may be defined in an operational sense by first performing 
all computations with the ordinary QFT dual to a holographic RG flow, and then continuing the 
number of colors N and the momenta appropriately. In the large N limit the correlators are given 
as a series in 1/-/V 2 , and so the analytic continuation simply amounts to inserting minus signs. 
Furthermore, the continuation in momenta is such that the effective coupling g^ s = gy M N/q, 
where q is a momentum, remains real. It would be very interesting to understand such 'pseudo'- 
QFTs from first principles. This would allow for a non-perturbative definition of the theory and 
would presumably elucidate many puzzling features of quantum gravity such as, for example, 
the entropy of de Sitter spacetime. Note also that the dual QFT provides a complete smooth 
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description of the system, including that of the initial singularity that is (generically) present 
in the background FRW solution. Our focus here was to extract the late-time behavior of the 
system, but it would be very interesting to understand the implications for singularity resolution. 

For the present, however, the main question is whether or not holographic models exist that 
are compatible with current observations but have a distinct phenomenology from standard 
inflation. We have shown that the answer to this question is affirmative. Initiating a 
holographic phenomenological approach, we found that it is straightforward to satisfy the 
current observational constraints using simple QFT models containing only a few parameters. 
The near scale invariance of the cosmological power spectra follows immediately from simple 
dimensional considerations. A number of parameters in the QFT model may then be estimated or 
constrained using the present observational data; once the remaining parameters have been fixed, 
all cosmological observables (including non-Gaussianities) then follow from direct computation. 
Note that these are complete models: there are no UV issues as these theories are super- 
renormalizable and furthermore the dimensionful coupling constant acts as an infrared regulator. 
Even without knowing the values of all the parameters in the QFT model, it is still perfectly 
possible to obtain concrete predictions, since not all cosmological observables depend on the full 
set of QFT parameters. Examples include the scale invariance of the power spectra at leading 
order, as well as the running of the spectral index discussed above. In general, we expect to 
obtain predictions that are qualitatively different from those of standard inflationary scenarios 
based on weakly coupled gravity. These expectations are borne out by the form of the running 
we found for the scalar spectral index. 

Clearly, the proposed phenomenological approach to holographic cosmology is worthy of 
further development. Over the next few years, forthcoming experiments (in particular the Planck 
satellite) promise to dramatically improve the observational constraints on many important 
cosmological parameters. It may well be that future observations confirm the predictions of 
holographic models of the type advocated here. The success of such an endeavor might then 
provide the first observational evidence for the holographic nature of our universe. 
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